Correlated phases of bosons in the flat lowest band of the dice lattice 
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We study correlated phases occurring in the flat lowest band of the dice lattice model at flux 
density one half. We discuss how to realize the dice lattice model, also referred to as the 75 lattice, 
in cold atomic gases. We construct the projection of the model to the lowest dice band, which yields 
a Hubbard-Hamiltonian with interaction-assisted hopping processes. We solve this model for bosons 
in two limits. In the limit of large density, we use Gross-Pitaevskii mean-field theory to reveal time- 
reversal symmetry breaking vortex lattice phases. At low density, we use exact diagonalization to 
identify three stable phases at fractional filling factors v of the lowest band, including a classical 
crystal at = 1/3, a supersolid state at = 1/2 and a Mott insulator at = 1. 



The dice lattice (TJ [2] , which is also referred to as the 
73 lattice [2, gives rise to an exceptional mechanism for 
localization: when subjected to a magnetic field with a 
flux density of one half flux quantum per plaquette, 
it realizes "Aharonov-Bohm cages" that perfectly con- 
fine particle motion by destructive interference around 
individual plaquettes [3]. The presence of these local- 
ized states translates into a macroscopic degeneracy of 
states or equivalently a spectrum with flat energy bands. 
A natural implementation of this model can be found 
in Josephson junction arrays [4 , where condensates of 
Cooper-pairs on each site are well characterized by an xy- 
model for the classical order parameter. The rich physics 
of these systems stems from a highly degenerate manifold 
of states at low energies [5 , described as vortex lattices. 
Dynamics in this manifold is slow [6 , and ordering at the 
lowest temperatures is determined by such subtle effects 
as magnetic interactions of currents [7| or anharmonic 
fluctuations [8 . 

Given the perfectly flat bands in the spectrum of the 
dice lattice, interactions can potentially lead to strongly 
correlated states in the regime of low particle density 
where number fluctuations are significant. This is remi- 
niscent of the physics of the fractional quantum Hall ef- 
fect [9 , which is predicted to exist also in systems of 
bosonic atoms in the continuum or in the Hofs- 

tadter bands of a square lattice pTHTS] . Similarly for 
fermions, there has been intense interest in realizing frac- 
tional quantum Hall states in general flat-band models 
with nonzero Chern numbers [14 . However, the dice lat- 
tice model at flux density of one half is fundamentally 
different from this Landau level physics, as it does not 
break time-reversal symmetry. In this model, the role 
of interactions has been described only for the two-body 
problem of fermions, where it leads to delocalized zero- 
energy states of pairs [I5j [16] . 

In this paper, we study the physics of strongly cor- 
related states in the flat lowest band of the dice lat- 
tice model at flux density = 1/2. For additional 
motivation of this study, we sketch candidates for ro- 
bust implementations of the dice lattice model in cold 
atoms. In constrast to other flat band models, as on 



the kagome lattice [17 , our proposal benefits from the 
magic of Aharonov-Bohm cages and thus yields a per- 
fectly flat lowest band separated by a large gap from 
higher bands, using only nearest-neighbor hopping. To 
analyse the problem of interacting particles in this flat 
band, we introduce an effective model obtained by pro- 
jecting the Hamiltonian to this band, which is character- 
ized by one dimensionless interaction parameter. In the 
limit of large density, we solve this model using Gross- 
Pitaevskii mean-field theory and find time-reversal break- 
ing ground states that realize the same phase patterns as 
in the ground states of the classical x?/- model [5 , but 
which have additional density modulations. The pro- 
jected model also enables efficient numerical studies of 
the many body physics on the dice lattice based on ex- 
act diagonalization. We undertake a numerical study in 
the regime of low particle density that is of interest for 
cold atomic gases on lattices. We identify stable phases 
of the model as a function of density and find evidence 
for several phases including classical crystalline states, a 
supersolid state with triangular crystalline order and a 
Mott insulator. 

An optical dice lattice can formed by three mutually 
phase coherent pairs of counterpropagating laser beams 
at relative angles of 27r/3 [18 . Here, we propose instead 
to use a set of three mutually incoherent pairs in the 
same geometrical arrangement. This enables a particu- 
larly convenient scheme to be realized for Yb by using 
an "antimagic" wavelength, such that two internal states 
(^5'o and ^Pq) are trapped at the points of maximum or 
minimum laser intensity these points again form a 
dice lattice, with the structure shown in Fig.[l] To induce 
dynamics in this lattice, one employs laser- assisted hop- 
ping that simultaneously imprints phases onto the hop- 
ping matrix elements [19]. For the flux density, = 1/2, 
that yields the desired flat band structure on the dice lat- 
tice [3 , the magnetic unit cell contains six inequivalent 
atoms. Following the ideas of [19 , the required phases 
can be imprinted by optically assisted tunnelling involv- 
ing several coupling lasers and one additional superlat- 
tice laser to break inversion symmetry of the magnetic 
unit cell. We note that at the field strength = 1/2 
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FIG. 1. (Color online) Figure showing the structure of the 
dice lattice, highlighting a rectangular unit cell and negative 
hoppings on three bonds (double hatched) that realize the 
case of flux density = |. The blue contour highlights the 
extent of a maximally localized single-particle wavefunction 
in the flat lowest band. Localized states centered around all 
six-fold connected sites of the lattice (encircled) form an or- 
thonormal basis for this band, and lie on a triangular lattice. 



time-reversal symmetry is not broken, so it is possible to 
choose a gauge with real tunnelling amplitudes and with 
only three bonds in the unit cell having a negative sign, 
as indicated in Fig. [l] This allows for an alternative 
potential implementation, using an optical dice lattice 
for single species atoms [18 . Negative hopping can be 
achieved in principle by shaking lattice sites [20 ; the ap- 
plication in this case may be challenging, as it requires 
the three-fold connected sites in the magnetic unit cell 
to be independently shaken. Finally, the synthetic gauge 
field can also be simulated by the Coriolis force due to a 
rapid rotation of the system [211 ES] . 

Our main focus concerns the many-body physics of 
bosons on the dice lattice at this flux density, = 1/2, 
as described by the Bose-Hubbard model 

n = -t {dlcLje'^^^ + h.c.) (1) 
<i,/i> 

/iGA,V jG* 

Here, a^^"^ are the creation/annihilation operators of 
atoms on single lattice sites. We adopt Roman indices 
for six-fold connected 'hubs' (*) and use Greek indices 
to number threefold (A, V) connected sites of the dice 
lattice. These types of site can have distinct on-site re- 
pulsion /7*, /7a. We choose a real gauge with A^j = or 
TT, with negative bonds highlighted in Fig. [l] 

The single-particle spectrum of dl]) at = 1 /2 is given 
by three flat bands with E = — VoT, 0, ^/6t ^. We focus 
on the regime with nU^/i^/t <C 1, such that all dynamics 
occurs in the lowest band only. Wavefunctions localized 
on the hubs and extending to the neighboring threefold 
connected sites (see Fig. [T]) span an orthonormal basis 



{^j} for the lowest band p!6l [23] , These <j)j have am- 
plitude l/>/2 on the central site j, and l/>/T2 exp[iAj^] 
on the peripheral sites /i [24]. In our real gauge, are 
identical up to translation in all unit cells. (Both sublat- 
tices are also related by a magnetic translation involving 
a gauge transformation.) The projection onto the low- 
est band realizes a new effective problem on a triangu- 
lar lattice visualized in Fig. [2| which derives from the 
density- density interactions of the microscopic Hamil- 
tonian Q via Vijki = (m)</>/c(m)0/(m) + 

^* Y.n^^)^)(^)^k{q)(t)M), and that reads 



=71 rii{ni - 1) + 72 XI [^^^j + + cfcj 
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with c\ creating an atom in orbital ^j, and 



71 = + ^/7a, 72 = ii4^A, 73 

A 



actors & 



qr 



St 



48 

= ex.pi[A 



- C-Cj, 

^/7a. The pref- 



288 



lis 



A^t] derive from 
the gauge fields around a shared threefold connected site 
fi [see Fig. [ij^b-vi)]. In the real gauge, cr]^^ = aj^ re- 
duces to a sign as illustrated in Fig. [2|c). This effective 
model includes on-site and nearest-neighbor interactions 
as well as pair hoppings, and interaction-assisted hopping 
processes. Equivalent projected models can be written 
generally also for multicomponent bosons/fermions. 

Importantly, note that while the original model ([T]) is 
strictly local, the projected Hamiltonian involves nonzero 
range interactions between nearest-neighbor sites. In 
contrast to the flat-band physics of the kagome lattice 
[17], longer range interactions vanish exactly. The ratio 
of on-site terms to nearest-neighbor terms can be tuned, 
since u = U^/U^ remains as a free parameter of the 
model (while the tunnelling t drops out of the problem). 
Choosing u ^ 1 while maintaining nU^^i^ <C t defines 
the hardcore limit of the projected model, eliminating 
configurations involving double occupancy. 

We now inquire into the nature of the many-body 
ground states of spinless bosons as a function of den- 
sity. At densities n ^ 1, but remaining in the regime 
n/7*^A <^ ^ governed by the projected Hamiltonian ([2|, 
we expect that correlations can be neglected. We an- 
alyze this regime using a Gross-Pitaevskii mean-field 
equation deriving from the matrix elements Vijki in 
the lowest band. This ansatz is based on a conden- 
sate wavefunction which is a coherent state of the form 
1^) = expE_^. a_^-c]]|0), such that c_^-|^) = In- 
troducing a chemical potential /i, the energy is {H) = 
\ T^ijki Vijkia*a*akai - /j.^. Using steepest de- 

scent numerical minimization of the energy, and expand- 
ing the resulting states on the full dice lattice, we find 
that the phase-patterns minimizing the energy are iden- 
tical to those described by Korshunov for the xy-model 
[5 . These ground states break time-reversal symmetry 
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FIG. 2. (Color online) (a) The projective dynamics within 
the low energy band of the dice lattice at flux = 1/2 re- 
alizes an effective model on a triangular lattice with an en- 
larged magnetic unit cell of two distinct triangular sites A and 
B. (b) The density- density interactions of the dice Hamilto- 
nian M give rise to the five distinct processes i)-v), including 
an on-site (i) and nearest-neighbor (ii) interactions, coherent 
hopping of pairs of particles onto the same (iii) and two dis- 
tinct neighboring sites (iv), and a stimulated hopping process 
(v). Hoppings are shown as arrows, density interactions as 
dotted lines. Processes (iv) and (v) have a nontrivial depen- 
dency on the gauge from (vi), see main text, c) Signs for the 
processes (iv,v) in the real gauge of Fig. [l] these processes in- 
volve sites on a triangular plaquette, of which there are four 
types classified by the participating sublattice indices [empty 
(A) / full (B) circles]; in addition one of the three sites carries 
two creation/annihilation operators (marked by a hexagon). 



and feature clusters of three plaquettes sharing the same 
vorticity or phase winding ±7r. Unlike for the x?/- model, 
the density for these mean- field states of the projected 
model ([2| is not homogeneous. Instead, we find = 2nA 
for 1^ = 1. Quantum fluctuations are likely to break the 
degeneracy of these different vortex patterns, as occurs 
for thermal fluctuations in Josephson junction arrays [8 . 
However, this effect is beyond Gross-Pitaevskii theory, 
which does not capture the role of quantum fluctuations. 

To address the question of how fluctuations resolve the 
frustration seen in the mean-field solutions, we study 
the regime of low particle density that takes the sys- 
tem into the strongly correlated limit. We express n 
in terms of the density of states in the lowest band n^, 
defining the band filling factor u = n/ug = 3n. A re- 
markable feature of the Hamiltonian ([2| is that all hop- 




FIG. 3. (Color online) (a) Energy per particle as a function of 
the filling factor = 3n of the low- lying band for the projected 
Hamiltonian ([2]) at = 1/2 with hardcore interactions. Fi- 
nite size effects between different lattice geometries and the 
effect of twisted boundary conditions (band-width shown as 
error bars for 3x3 and 3x4 geometries) are small. Inset: En- 
ergy measured with respect to the tangent at points = 1/2 
and zy = 1 (dashed line in main panel), showing that = 1/2 
is the only stable phase between the crystal at zy = 1/3 and 
Mott state at z/ = 1 for u ^ 1. (b) Second difference of the en- 
ergy per particle. Maxima indicate potential incompressible 
phases, including the states at u = 1/3, 1/2 and 1 as well as 
an additional weaker candidate at z/ = 2/3 that is overridden 
by phase separation. 



pings are mediated only by the presence of neighboring 
particles, while single particles remain stationary. Con- 
sequently, the many-body spectrum at low filling factor 
u < = 1/3 has a highly degenerate zero-energy ground 
state. The two-body problem has been studied for two- 
species fermions with repulsive contact interactions [16], 
and it features delocalized zero-energy spin-singlet states 
in addition to the trivial zero-energy states with parti- 
cles placed at a distance of more than one lattice vector. 
The two-boson wavefunctions have the same structure, 
due to their identical spatial symmetry. Counting ^ = 
states in the spectra of finite size systems confirms this. 
For filling factors u < Uc = ^/^ the system has infinite 
compressibility. Precisely at the model yields an in- 
compressible ground state that is a classical crystal of 
triangular symmetry with unit vectors Ui = ffi -\- f]2 and 
U2 = 2f]i — f]2 (see Fig. [l] for the definition of 7/1,2)- The 
three degenerate crystal ground states can be written as 
l^c) = Yln,m^^[^t + nui + mu2]\0), whcrc c^[r\ fills the 
orbital centered at r, and ft ^ {0,7/1,7/2} are transla- 
tions with respect to the origin. Compression to densities 
u > Uc = 1/3 costs a nonzero energy. Such incompress- 
ible states at noninteger filling can exist only due to the 
nonlocal interactions of the effective model. 

To study the behaviour of the model at densities 
u > 1/3, we calculate its spectrum and eigenstates using 
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FIG. 4. (Color online) Data for N — 8 atoms at ly — 1/2 on a lattice of 4x 2 unit cells, (a) Two-point correlation function {hrfio) 
of the groundstate |0) for — (0,7r). The reference site '0' is visually highlighted. Axes carry units a — \vi\^ h — ^3/4 |. 
(b) Density {fir) of the symmetry-broken crystal state \S) (see main text), (c) Fluctuations (n^) — (rir)^ (blue circles), and 
condensate wavefunction vq of \S) (red arrows), (d) Scaling of the condensate fraction \q/N against A^~^, for superposition 
states \S) and pinned states with pinning centers of strength Vp on crystal sites. 



exact numerical diagonalization, focusing on the hard- 
core limit u ^ 1. We consider a finite size system con- 
sisting of Ni) bosons in x Ly rectangular unit cells 
of the structure shown in Fig. |2|a). We consider peri- 
odic boundary conditions allowing arbitrary flux (Ox^Oy) 
passing through the two cycles of the ensuing torus. 

We first ask at which densities there can appear phases 
that are stable to phase separation. This requires positive 
compressibility, such that A'^E{N) = E{N ^1) ^ E{N - 
1) - 2E{N) > 0. The dependency E{iy) for < < I 
is shown in Fig. ^di) and A'^E in Fig. jsj^b). Besides 
the crystal described above at Uc = 1/3, three further 
densities u = 1/2, u = 2/3, and u = 1 emerge as natural 
candidates for stable phases, signalled by a sharp peak in 
A'^E. Our results for the energy per particle in Fig.jsja) 
indicate that the state of intermediate density u = 2/3 
is merely a local minimum and is unstable to phase sep- 
aration into regions of density u = 1/2 and 1 at long 
timescales, as demonstrated in the inset by a Maxwell 
construction for the energy of a phase separated system 
composed of the adjacent stable phases. The physics of 
the V = 1 state is readily identified as that of a Mott 
insulator. 



The presence of the additional stable state at u = 1/2 
is an intriguing feature of the projected model Our 
exact diagonalization studies at this filling, for systems 
up to 16 particles, reveal a considerable dependency on 
the boundary conditions 0. The eigenvalues disperse, and 
bands cross, so there is no protected groundstate mani- 
fold. Nonetheless, at the energy minimum, [occurring at 
e = (0,7r), (7r,7r), and (0,0) for 4 x 2, 6 x 2, and 4 x 2 
lattices, respectively] the lowest four eigenvalues are sep- 
arated from higher excited states by an amount larger 
than their respective splitting. The correlation function 
of the lowest eigenstate at this point [Fig. [4]^a)] reveals a 
triangular crystalline order with lattice vectors {2iTi, 2v2} 
that suggests a superlattice cell composed of 4 sites, con- 
sistent with the relative isolation of the four low-lying 
eigenstates. 

Taking a superposition of the four low-lying (and trans- 
lationally invariant) exact eigenstates, one can explicitly 
construct symmetry-broken states [25 . The picture that 
emerges from these states is that of a supersolid, with 
half the particles forming a crystal while the remaining 
particles Bose condense in a state that fills the channels 
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between the former. With hindsight, the presence of a 
condensate fraction in the system explains the stiffness of 
the state to a twist in the periodic boundary conditions 
6 observed in the exact spectra. 

The density matrix p for a simple crystalline state of 
all Ni) bosons would be characterized by Ni) eigenvalues 
of order one. For the superposition states, we typically 
find one eigenvalue, Aq, that is larger than one and A^^,/2 
eigenvalues of order one, signalling the presence of both 
superfluid and crystalline components. Generalizing the 
construction of symmetry-broken condensate states in 
Ref. [25 , we generate these supersolid states as the su- 
perpositions of the four lowest eigenstates in the exact 
spectrum by maximizing = Y1\Zq ^^e sum of the 
first k = N/2 + 1 eigenvalues of p. In Fig.jljb), we show 
the density profile for the symmetry-broken state \S) ob- 
tained by optimizing P(n/2-\-i)^ with obvious crystalline 
order. In this state, only half the density is concentrated 
on superlattice sites, while the remainder forms a con- 
densed background liquid. On the crystal sites, particle 
number fluctuations are negligible, and there is a very 
small amplitude for the eigenvector of the largest den- 
sity matrix eigenvalue on these sites [Fig. [4]^c)]. This 
eigenvector has support on the background sites, where 
simultaneously fluctuations are strong. These features 
are characteristic of a supersolid fTT . We analyze Aq 
for the superposition states \S) constructed above, and 
alternatively for eigenstates obtained using trapping po- 
tentials Vp < on the superlattice sites. Indeed, our 
numerics indicate the condensate fraction extrapolates 
to a nonzero value 0.06 in the thermodynamic limit 
[Fig. |4jd)]. Unlike the vortex lattices at high density, we 
find that the condensed fraction of this supersolid does 
not break time-reversal symmetry. Experimentally, such 
crystalline order could be observed clearly in in situ im- 
ages of the system; expansion images should reveal evi- 
dence for a nonzero condensate fraction, and the density 
ordering will appear in the noise correlations [26^ . 

In summary, we have proposed the dice-lattice model 
as an attractive experimental system to explore the 
many-body physics of time-reversal symmetric flatbands. 
We have derived the effective Hamiltonian for interacting 
particles in such a band. Our results give the first insights 
into its rich phase diagram for bosons, driven entirely by 
interaction-assisted hopping processes which emerge as a 
generic feature of projected flat-band Hamiltonians. 
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